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ABSTRACT 

Context. Little is known about the planetary systems around single white dwarfs although there is strong evidence that they do exist. 
Aims. We performed a pilot study with the extreme-AO system on the Spectro-Polarimetric High-contrast Exoplanet REsearch 
(SPHERE) on the Very Large Telescopes (VLT) to look for giant planets around a young white dwarf, GD 50. 

Methods. We were awarded science verification time on the new ESO instrument SPHERE. Observations were made with the In- 
fraRed Dual-band Imager and Spectrograph in classical imaging mode in H band. 

Results. Despite the faintness of the target (14.2 mag in R band), the AO loop was closed and a strehl of 37% was reached in H band. 
No objects were detected around GD 50. We achieved a 5-sigma contrast of 6.2, 8.0 and 8.25 mags at 0'.'2, 0'.'4 and 0'.'6 and beyond, 
respectively. We exclude any substellar objects more massive than 4.0 Mj at 6.2 AU, 2.9 Mj at 12.4 AU and 2.8 Mj at 18.6 AU and 
beyond. This rivals the previous upper limit set by Spitzer. We further show that SPHERE is the most promising instrument available 
to search for close-in substellar objects around nearby white dwarfs. 
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1. Introduction 


Radial velocity and transit surveys demonstrate that extrasolar 
planets are prevalent around main-sequence stars in the Milky 
Way (Gumming et al. 2008). However, little is known about the 
evolution of planetary systems beyond the main sequence life 
time of their host stars. Planets have been discovered around 
(sub)giant stars (Johnson et al. 2011; Gettel et al. 2012) as 
well as horizontal branch stars (Silvotti et al. 2007; Setiawan 
et al. 2010; Charpinet et al. 2011). There is only one planetary 
mass object detected with direct imaging around a white dwarf 
(WD 0806-661; Luhman et al. 2011) although theoretical calcu¬ 
lations show that many more planets can survive the red giant 
stage (Burleigh et al. 2002; Jura 2008; Mustill & Villaver 2012; 
Nordhaus & Spiegel 2013). 

From an observational perspective, the strongest evidence 
for planets around single white dwarfs comes from the study of 
heavy-element-enriched white dwarfs. About 25-50% of white 
dwarfs have heavy elements in their atmospheres, in addition 
to the primordial hydrogen or helium (Zuckerman et al. 2010; 
Koester et al. 2014). There must be a constant supply of material 
onto the white dwarf to balance the rapid settling of heavy ele¬ 
ments out of the photosphere (Koester 2009). In addition, some 
of these heavy-element-enriched white dwarfs also display in¬ 
frared excess from a circumstellar dust disk (Jura et al. 2007; 
Xu & Jura 2012), or more rarely a gaseous disk (Gansicke et al. 


* Based on observations made with European Southern Observatory 
(ESO) telescopes at the La Silla Paranal Observatory under program 
60.A-9373(A). 


2006). The standard model is that these were formed by aster¬ 
oids that were perturbed by the surviving planet(s) in the sys¬ 
tem onto orbits that crossed the tidal radius of the white dwarf, 
were disrupted, and formed a dust disk; eventually, these materi¬ 
als were accreted onto the white dwarf, enriching its atmosphere 
(Jura 2003; Jura & Young 2014). This asteroid accretion model is 
further supported by comparing the abundances of the accreting 
material with solar system objects - they all resemble terrestrial 
material rather than the cosmic abundance (Klein et al. 2010; 
Jura et al. 2012; Gansicke et al. 2012; Xu et al. 2014). 

Dynamical simulations show that there are different path¬ 
ways that can lead to the asteroid accretion model. The config¬ 
uration varies between having a comet belt with two planets - 
an inner planet and a Neptune-sized planet (Bonsor et al. 2011), 
an asteroid belt with a Jupiter-sized planet (Debes et al. 2012) 
and an asteroid belt with an eccentric sub-Jupitzer-sized planet 
(Frewen & Hansen 2014). At least one planet must be present to 
explain all the observed features of the system. 

However, classical planet-hunting methods do not work well 
for white dwarfs. The red giant star has cleared planets within a 
few AUs, making the transit probability very low. White dwarfs 
have far too few spectral lines to make accurate radial velocity 
measurements. Here, we summarize a few previous searches for 
planets arounds single white dwarfs. 

- Direct imaging. It has long been speculated that direct imag¬ 
ing should be the most promising method to detect plan¬ 
ets around white dwarfs due to; (i) the better contrast be¬ 
tween the white dwarf and the potential planet compared to 
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main sequence stars; (ii) the orbital expansion of the planet 
during the red giant stage (Burleigh et al. 2002; Gould & 
Kilic 2008). The only planetary mass object around a single 
white dwarf was discovered via direct imaging with Spitzer; 
WD 0806-66IB has a mass of 7 Mj at an orbital separation of 
2500 AU (Luhman et al. 2011). The most systematic search 
is the DODO (Degenerate Objects around Degenerate Ob¬ 
jects) survey, which observed 29 white dwarfs with Gem- 
ini/NIRI and VLT/ISAAC; their average upper limit is ~ 8Mj 
beyond 35 AU (Hogan et al. 2011). There are also a few other 
direct imaging studies using HST/NICMOS (Debes et al. 
2005a; Zinnecker & Kitsionas 2008). 

- Infrared excess. Because a white dwarf’s luminosity peaks 
in the ultraviolet or optical, a cool companion can be easily 
detected as an infrared excess. Those searches started in the 
late 80s and so far, only brown dwarf companions have been 
discovered (Zuckerman & Becklin 1987; Maxted et al. 2006; 
Farihi et al. 2008; Kilic et al. 2009). 

- Planetary transit. Faedi et al. (2011) studied 194 white 
dwarfs in the WASP survey and put an upper limit of 10% 
for Jupiter-sized brown dwarfs and giant planets with orbital 
periods less than 0.1-0.2 d. Fulton et al. (2014) searched for 
eclipses around ~ 1700 white dwarfs in the Pan-STARRl 
field and found zero candidates, giving an upper limit of 
0.5% for Jupiter-sized-planet orbiting just outside of the 
Roche limit. 

- White dwarf pulsation. Some white dwarfs have extremely 
stable pulsation modes and the presence of a planet can cause 
a periodic change in the arrival time of the pulsation. Mul- 
lally et al. (2008) led a study on 15 such stable pulsators 
and identified one candidate - a 2Mj mass planet in a 4.5 yr 
orbit around GD 66. The planet hypothesis was dismissed 
later due to some other peculiar behavior in the pulsation 
of GD 66 (Mullally et al. 2009; Farihi et al. 2012; Hermes 
2013). The upper limit from pulsation studies is typically a 
few Jupiter masses within 5 AU of the white dwarf (Mullally 
et al. 2008). 

In this paper, we present a pilot search for giant planets 
around a young white dwarf GD 50. The data were taken in a 
science verification (sv) program with the new extreme-AO in¬ 
strument SPHERE (Beuzit et al. 2008) at ESO’s VET. SPHERE 
includes the SAXO extreme adaptive optics system (Fusco et al. 
2014; Petit et al. 2014), with a deformable mirror of 41 x 41 
actuators, and 4 control loops (tip-tilt, high-orders, near-infrared 
differential tip-tilt and pupil stabilization). 

2. Observation and Data Reduction 

2.1. Target Selection 

GD 50 (WD 0346-011) is a young nearby ultramassive hydro¬ 
gen atmosphere white dwarf (my = 14.0, mR=14.2, mn = 14.9). 
Its formation scenario has puzzled many astronomers. Dobbie 
et al. (2006) performed a complete astrometric and spectroscopic 
analysis and concluded that it has evolved as a single star and is 
an escaping white dwarf from the Pleiades open cluster. By fit¬ 
ting the photometric and spectroscopic data, they derived the ef¬ 
fective temperature, surface gravity and distance of GD 50 to be 
T = 41550 K, log g = 9.15 and d = 31 pc. This corresponds to a 
white dwarf with a mass of 1.264 Mq and a total age of 125 Myr, 
consistent with the age of Pleiades. According to the initial-final 
mass relation in Williams et al. (2009), it had an initial mass of 
7 Mq, making its progenitor a B-type star. 


GD 50 was chosen as our target for the SPHERE sv program 
due to its youth and close proximity. At 125 Myr, potential plan¬ 
ets can still be quite hot and thus detectable. In addition, GD 50 
has been observed with Spitzer/IRAC and its 3 cr upper limit 
of ~ 4 Mj ' is among the best constraints from Spitzer searches 
(Earihi et al. 2008). A side goal of the project is to compare the 
performance of SPHERE with other instruments for detecting 
substellar objects around white dwarfs. 

2.2. Observation 

GD 50 was first observed on December 5th, 2014. Due to poor 
observing conditions, the data were not used for the analysis. 
The observation was repeated on Eebruary 1st, 2015 with seeing 
around I'.'O, a moon distance of 46° (89% illumination) and a 
coherence time of ~ 5 ms. The InfraRed Dual-band Imager and 
Spectrograph (IRDIS) (Dohlen et al. 2008) was used in classi¬ 
cal imaging (Cl) mode in broadband H. The observations were 
made in pupil-stabilized mode, ideal for PSE subtraction, as the 
pupil and PSE stay fixed while the field rotates with respect to 
the detector. The integration time was 64 sec per exposure and 
the observing block consisted of a 4 x 4 dither pattern with 3 
exposures at each position, giving a total on-source time of 3072 
sec. 

2.3. Data Reduction 

The raw images were first corrected for Hatfield errors and bad 
pixels. The white dwarf was unsaturated in all the images. Since 
no coronagraph was used, we were able to align the star on all 
images with subpixel accuracy. Due to the faintness of the star, 
the stellar halo is weak and there are hardly any speckles. Thus 
there is no need for PSF-subtraction. All images were derotated 
to align North up and median combined. Eollowing Wahhaj et al. 
(2013), a filtering process was applied to the median-combined 
image to isolate spatial features of ~ 50 mas in size, the diffrac¬ 
tion limit in H band. This filter removes azimuthally extended 
features, e.g. the stellar halo, and improves the sensitivity to 
point sources for separations less than 073. The signal-to-noise 
map is shown in Fig. 1. 

Then we calculated the 5cr contrast curve, or the sensitivity 
limit relative to the primary. The peak intensity (star peak) at 
the stellar centroid was estimated by cubic interpolation. The 
pixel-to-pixel RMS was then measured in annuli of 5 pixels radii 
centered on the star for all separations (radii) of interest. The 5cr 
contrast was computed as (star peak)/(RMS x 5.0), where RMS 
is a function of separation from the star. No sources except the 
white dwarf were detected and the contrast curve is shown in 

Fig. 2. 

To check the robustness of our sensitivity limit, we cut out 
a circular aperture of radius 4 pixels around the star from each 
basic reduced image and re-inserted it at the 5cr contrast limit at 
several separations and position angles in the individual images. 
All the images were derotated and median-combined again. As 
shown in Fig. 1 , all the inserted companions are recovered in 
the median combined image. We have achieved 5cr contrasts of 
6.2, 7.8, 8, 8.1 and 8.25 mags at 072, 073, 074, 075 and 076, 
respectively. In the background limited region, the 5cr limit cor- 

* Farihi et al. (2008) reported a 3 cr upper limit of 3 Mj assuming a 
total age of 100 Myr for GD 50. To make the results directly comparable 
with our SPHERE observation, we adopted the revised age of 125 Myr 
and the updated BT-Settl model grid (Allard 2014). Our updated Spitzer 
upper limit is 4.2 Mj for GD 50. 
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responds to an apparent magnitude in H band of mn = 23.15 
(absolute magnitude Mh -20.7). 


• X 



1 arcsec 


Fig. 1: Signal-to-noise map (i.e., reduced image divided by a 
noise map derived from it, see Section 2.3). The map has been 
smoothed using a Gaussian of FWHM 1.5 pixels. No astronom¬ 
ical point-like features were detected. In addition, we recovered 
artificial sources inserted into the raw data at the 5cr contrast 
level at separations ranging from 100 mas to 3 arcsec with in¬ 
crements of 100 mas (X shaped pattern of sources). The held of 
view displayed here is 6'.'25 by 6"25. Counts are displayed in lin¬ 
ear scale and inverted (dark corresponds to high counts). North 
is up. The central white dwarf is not visible in this map because 
of image processing. 



Fig. 2: 5cr contrast curve for GD 50. A 5(T contrast of 6.2, 7.8, 
8.0, 8.1 and 8.25 mags has been reached at 200 mas, 300 mas, 
400 mas, 500 mas and 600 mas, respectively. 


aration as the orbital distance from the white dwarf and did not 
correct for the inclination effect. As shown in Fig. 3, we put an 
upper limit to the planetary mass of 4.0 Mj, 2.9 Mj and 2.8 Mj 
at 6.2 AU, 12.4 AU and at 18.6 AU and beyond, respectively. 

Our upper limit corresponds to a limit of 2.8 Mj beyond 
3.4 AU during the main-sequence stage of the star. Depending 
on the metallicity, a TMq star could reached a maximum radius 
about 6.5 AU during the post main-sequence evolution (Veras 
et al. 2013), so any planets within 6.5 AU would possibly have 
been destroyed. Therefore, it is unlikely that there are any giant 
planets more massive than 2.8 Mj around GD 50 now. For the 
surviving planets, their orbits will expand by a factor of about 
7/1.264, which will be at least 36 AU from GD 50. However, 
tidal interaction between the planet and the extended envelope 
of the AGB star can counteract the orbital expansion (Mustill 
& Villaver 2012). Furthermore, if multiple planets were present, 
mutual interactions could scatter planets inward during the white 
dwarf stage (Veras & Gansicke 2015). As a result, some planets 
might end up having a much smaller semi-major axis. In addi¬ 
tion, no heavy elements have been detected in the atmosphere 
of GD 50; though optical observations are not very sensitive to 
heavy elements at high stellar temperatures (Koester et al. 2009). 
GD 50 has a temperature of 41550 K and the upper limits are not 
very constraining. 

In Fig. 3, we compare the upper limit on GD 50 with all 
previous imaging searches for substellar objects around white 
dwarfs. Our upper limit on GD 50 is the most stringent, both in 
terms of mass and separation. For example, HST/NICMOS has 
an inner working angle of 079 with the coronagraph (Debes et al. 
2005a), which for GD 50 corresponds to 28 AU in Fig. 3. 



Fig. 3: Upper limits from previous imaging surveys (DODO sur¬ 
vey from Hogan et al. 2011, HST/NICMOS studies from Debes 
et al. 2005a) for substellar objects around white dwarfs. The only 
positive detection so far is WD 0806-66IB, but it is not shown 
because the separation is 2500 AU (Luhman et al. 2011). We 
also include the upper limit on GD 50 from Spitzer/IRAC stud¬ 
ies (Farihi et al. 2008). GD 50 is the most distant white dwarf in 
the figure, yet our limit on GD 50 is still the most stringent. 


3. Discussion 

3.1. Upper limits on planetary mass companions of GD 50 

To convert our detection limit to planetary mass, we used the 
BT-Settl models from Allard (2014). We took the projected sep- 


3.2. Future observational prospects 

GD 50 is not representative of nearby white dwarfs - it is partic¬ 
ularly young. Now we investigate the performance of SPHERE 
compared to other instruments for a generic white dwarf within 
20 pc, with an absolute magnitude My = 12.0, Mh = 12.0 and 
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a total age of 2 Gyr (Holberg et al. 2008). Since there are no 
BT-Settl models available for low-mass planets at 2 Gyr, AMES- 
COND models are used instead (Allard et al. 2001). As shown 
in the left panel of Fig. 4, for nearby white dwarfs (within ~ 30 
pc), Spitzer photometric observation is limited by the calibration 
uncertainty of approximately 15% (3cr); as a consequence, the 
Spitzer upper limit to planet mass remains the same regardless 
of the white dwarf distance. This is not the case for SPHERE 
because a potential planet would be spatially disentangled from 
the star outside an angular separation of ~ 0"2, about 2 AU at 
10 pc. Resolved imaging with Spitzer can be more sensitive to 
lower mass objects, but at a much larger orbital separation (e.g. 
WD 0806-661B in Luhman et al. 2011). Spitzer observations 
become background limited for more distant white dwarfs (be¬ 
yond 30 pc) and all those white dwarfs are typically too faint 
for SPHERE. For nearby white dwarfs, resolved imaging with 
SPHERE is more sensitive to close-in low mass objects than 
photometric observations. There are ~ 50 white dwarfs brighter 
than 14.5 mag in R band that are observable with SPHERE (Hol¬ 
berg et al. 2008). They are prime targets for future searches for 
planetary systems with SPHERE. 

In the right panel of Fig. 4, we show the performance of 
SPHERE compared to other direct imaging instruments^. NIC- 
MOS and NACO have a much bigger inner working angle. Note 
that the difference between NACO, HST, and SPHERE in the 
background limited regime is minor. The real improvement of 
SPHERE over other instruments is in the region less than I'.'O. 
For a typical white dwarf within 20 pc, this corresponds to within 
20 AU from the star. This region is particularly interesting, be¬ 
cause it can constrain models of the survival of planetary systems 
during the AGB phase. 

Most white dwarfs have a progenitor mass between 2-5 Mq. 
Searching for planets around white dwarfs directly probes this 
intermediate-mass region, which is very difficult with radial ve¬ 
locity and transit methods (Galland et al. 2005). Recently, Barber 
et al. (2014) discovered a 1.04 Mq white dwarf (progenitor mass 
5.4 Mq) with a dust disk from tidally disrupted asteroid(s) - the 
hrst confirmation that massive stars can have planetary systems. 


4. Summary 

We present a pilot study with the extreme-AO instrument 
SPHERE to search for giant planets around a young white dwarf, 
GD 50. We exclude any substellar objects above 4.0 Mj, 2.9 Mj 
and 2.8 Mj at 6.2 AU (0'.'2), 12.4 AU (0'.'4) and at 18.6 AU (0'.'6) 
and beyond, respectively - among the best constraints compared 
to all previous studies. We further show that SPHERE is the 
best instrument available to search for close-in substellar objects 
around nearby white dwarfs. 
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Fig. 4; Illustration of the performance of SPHERE compared to other instruments for a typical nearby white dwarf (age 2 Gyr, My 
= Mh = 12). The green regions represent the parameter space for which SPHERE enables significant improvement over the other 
instruments. Left: Sensitivity vs. distance of the white dwarf from the Sun. The quantity p represents the angular separation that the 
instrument is sensitive to. The performance of SPHERE is extrapolated from our observation of GD 50. Because GD 50 is at the 
faintness limit, the actual performance on a brighter target would be better. The vertical cut-off in the SPHERE curve near 30 pc is 
artificial and assumes that a star fainter than mv=14.5 cannot be observed due to the limiting magnitude of the AO. Eor nearby bright 
white dwarfs, resolved imaging is more sensitive to low mass objects than photometric observations. Right: Contrast curves for the 
same WD at a distance of lOpc. The NACO curve is estimated from archival data (program ID; 079.D-0561(A), PI; Radiszcz) 
with a similar integration time as SPHERE (3000 sec). The HST/NICMOS limit was extrapolated from Debes et al. (2005b) with 
a typical exposure time of 20 min. The key point is that SPHERE has the unique capability to detect high-mass planets at 2-10 AU 
from the white dwarf, which is not possible with any other instrument. 
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